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ABSTRACT 

The XA region of the Cygnus Loop is a bright knot of X-ray emission on the 
eastern edge of the supernova remnant. The emission resuhs from the interaction 
of the supernova blast wave with density enhancements at the edge of a precursor 
formed cavity. However, this interaction is complex given the irregular morphol- 
ogy of the cavity wall. To study the nature and origin of the X-ray emission we 
use high spatial resolution images from Chandra. We extract spectra from these 
images to analyze the physical conditions of the plasma. Our goal is to probe the 
density of various regions to form a picture of the cavity wall and characterize 
the interaction between this supernova and the local interstellar medium. We 
find that a series of regions along the edge of the X-ray emission appears to trace 
out the location of the cavity wall. The best fit plasma models result in two 
temperature component equilibrium models for each region. The low tempera- 
ture components have densities that are an order of magnitude higher than the 
high temperature components. The high density plasma may exist in the cavity 
wall where it equilibrates rapidly and cools efficiently. The low density plasma 
is interior to the enhancement and heated further by a reverse shock from the 
wall. Calculations of shock velocities and timescales since shock heating are con- 
sistent with this interpretation. Furthermore, we find a bright knot of emission 
indicative of a discrete interaction of the blast wave with a high density cloud 
in the cavity wall with a size scale ~0.1 pc. Aside from this, other extractions 
made interior to the X-ray edge are confused by line of sight projection of var- 
ious components. Some of these regions show evidence of detecting the cavity 
wall but their location makes the interpretation difficult. In general, the softer 
plasmas are well fit at temperatures (kT) ~0.11 keV, with harder plasmas at 
temperatures of (kT) ~0.27 keV. All regions displayed consistent metal deple- 
tions most notably in N, O, and Ne at an average of 0.54, 0.55, and 0.36 times 
solar, respectively. 
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Subject headings: supernova remnants — ISM: individual (Cygnus Loop) — X- 
rays: individual (Cygnus Loop) — plasmas — shock waves 



1. Introduction 



Supernova remnants are a dominant factor in the dynamics of a galaxy given that their 
ubiquity and size influence all phases of the interstellar medium (ISM). An investigation 
of these objects is crucial to understand this interaction and the mechanisms involved with 
transferring energy and matter to the ISM. The Cygnus Loop provides an excellent laboratory 
for the study of a galactic supernova remnant as it is large, ~ 3° x 3° ( Levenson et al.|[l997 ), 



close by, 540 pc (Blair et al. 2005), and relatively unabsorbed {E{B — V) = 0.08; Fesen et 



al. (1982)). Levenson et al. (1997, 1998) performed a detailed study of the optical and X-ray 



emission, and found that the morphology is not indicative of a typical theoretical blast wave 



in a uniform medium (Sedov 1993 Spitzer 1998). Instead, the precursor wind carved out a 



nearly spherical cavity made of clumps at high density separated by lower density gas, thus 
providing the ingredients for a cavity explosion ( McCray fc Snow||1979 ). Over the history of 
this remnant the blast wave has propagated through a low density plasma and has relatively 
recently interacted with the cavity wall leading to limb-brightened emission. The density 
enhancements decelerate the blast wave and lead to the X-ray emission. 

X-ray studies of specific blast wave-cloud interactions give mixed results. Several studies 
suggest low temperature (fcT ~0. 1-0.2 keV) equilibrium conditions within these dense clouds 



with higher temperatures toward the interior (Levenson et al. 2002 Levenson & Graham 



2005 Leahy 2004). Alternatively, Katsuda et al. (2008) and Miyata et al. (2007) find nearly 



all of the best fit temperatures above kT = 0.2 keV with nonequilibrium models. However, 
the general trend is a multiple temperature component model that is consistent with a lower 
temperature exterior and higher temperature interior, which may be additionally heated by 
reverse shocks arising from the cloud interactions. 

The XA region is a complicated, convoluted region of shocks located on the eastern 



side of the Cygnus Loop, on the south end of NGC 6992 (Hester & Cox 1986). Previous 



studies have looked at this region over a range of wavelengths (Sankrit et al. 2007 Danforth 



et al. 2001 Miyata & Tsunemi 2001). Sankrit et al. (2007) find complicated ionization 



structures within UV shocks. Danforth et al. (2001) also perform UV emission line studies 



and formulate a picture of the XA region where a finger of dense gas is protruding into 



this region. In addition, Miyata & Tsunemi (2001) observed this region with the ASCA and 
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ROSAT observatories. Similar to the above X-ray studies, they find inhomogeneous structure 
and explain their results with a refiection-shock model resulting from the interaction of the 



blast wave with a single cavity wall cloud. More recently, Zhou et al. (2010) analyzed higher 
resolution XMM-Newton data. They find a low temperature exterior {kT ~0. 07-0. 15 keV) 
with depleted abundances and high temperature interior {kT ~0. 24-0. 46 keV). They argue 
that the region is dense and clumpy, but identify no small scale features. 

The goal of this analysis is to exploit the high spatial resolution of Chandra to elucidate 
structure in the XA region and obtain more detail on the interaction of the blast wave with 
the precursor formed cavity wall. Combined with the ability to extract spectra, it will be 
possible to probe the XA region for additional complexities and attempt to model the cavity 
wall morphology and its physical state. 



2. Observations and Data Reduction 

The bright, shock rich region on the eastern edge of the Cygnus Loop, known as the 
XA region, was observed with the Advanced CCD Imaging Spectrometer on the Chandra 
X-ray Observatory for 75 ks. The observation was made on 2001 July 28 and archived under 
observation ID 1961. The region was centered on the S3 chip which is back-illuminated 
and therefore more sensitive to the large amount of soft X-rays being emitted from this 
region. This paper concentrates on the interaction of the supernova blast wave with the 
cavity boundary, which is dominated by soft X-rays. For that reason we have decided to 
concentrate our analysis on data contained within only the S3 chip giving a 8.3' x 8.3' 
field of view. The observation produced imaging with a spatial resolution of ~ 1" and 
spectral resolution E/AE~5-20 over an energy band of E~0.5-1.5 keV (fiux above 1.5 keV 
is negligible). 

The level 1 data were reprocessed to level 2 with standard processing procedures in the 



Chandra Interactive Analysis of Observations (CIAO, v4.1.2, Fruscione et al. (2006)) software 
package with current calibration data from the Chandra Calibration Database (CALDB, 
v.4.1.0). Good time intervals, charge transfer efficiency, and time dependent gain variations 
were accounted for. The data were background corrected using extractions from eastern 
areas of the chip that are devoid of emission. 

The reprocessed raw data is displayed in figure [T] It is evident that this region is 
extremely complex with many shocks existing down to the smallest scales. The emission 
appears to be concentrated in clumps or hot spots as opposed to filamentary shock struc- 
tures. To elucidate the spectral structure within this region, a false color image, created 
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by discriminating the data via energy, is shown in figure [2] The raw data are segmented 
into 0.3-0.5 keV (red), 0.5-0.8 keV (green), and 0.8-1.5 keV energy bands. The soft, red 
emission is dominant on the very eastern edge where the blast wave is initially encountering 
the cavity wall, although some red emission is seen toward the interior, especially on the 
south end. The harder, green emission is interior to this red emission as is the hardest, blue 
emission. To gain more insight into the distribution of these energy regions a ratio map 
of the soft emission (0.3-0.5 keV) to hard (0.8-1.5 keV) emission is displayed in figure [3] 
This shows the predominance of soft emission at the eastern edge where the blast wave is 
interacting with the cavity wall. Also, there is a large, bright, soft region centrally located 
in declination, along the front, that extends into the remnant considerably with relatively 
little high energy contribution. Finally, there is a line of soft emission that extends from the 
shock front westward into the interior of the remnant along a declination of ~ 31°01'. 



3. Spectral Extractions and Fitting 

We extracted data from specific regions in order to perform spectral analyses. These 
regions were chosen for their apparent physical qualities based on inspection of the false color 
image and the band ratio map. The combination of these images can be used to elucidate 
regions with distinct physical characteristics. Spectrally extracted regions are shown in 
figure |4] overlaid on both the false color image and the band ratio map. Spectra from 
these regions were created using the specextract tool in CIAO. This tool also creates specific 
ancillary response files and redistribution matrices for each region. 

The first region of interest is the soft X-ray dominated rim of the remnant. The extrac- 
tion regions are depicted by the series of rectangles following the rim and labeled B1-B5. 
Inspection of the color image shows that these regions are dominated by soft emission. In 
contrast to these B regions, the regions Al and A2 are dominated by harder emission with 
A2 located in a veritable void in the ratio map. The next region of interest is the centrally 
located soft X-ray enhancement, probed by extraction regions CI and C2. CI corresponds to 
a bright knot of emission in the false color image while C2 corresponds to a bright spot in the 
ratio map. This area of soft X-ray enhancement appears extremely complex with structures 
appearing more like clumpy hot spots rather than elongated shock filaments. Another region 
of interest lies along the westward soft X-ray extension along declination ~ 31°01'. Regions 
E1-E3 are used to investigate this emission. Finally, there is an extremely bright knot in 
the southeast located at ~ 31° of declination and right ascension of 20'^ 57mm 21'"^'^. The 
high count rates allow for several extractions over this region. A grid of extraction boxes is 
used to investigate the spectral properties of this bright emission. 
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To analyze the extraction regions we use the Xspec spectral fitting package vl2.5.1n 
(Arnaud 1996). Only data above 0.5 keV are fit due to the uncertainty of the low energy 



calibration. The data are also binned to > 30 counts per bin so that Gaussian statistics can be 
applied. If the X-ray emission arises from supernova blast wave heating of the ISM, then we 
expect to see cooling emission predominantly in emission lines. Depending on the density of 
the plasma and the time since it has been shocked, it will either been in coUisional ionization 
equilibrium (high density or large amount of time since heating) or on the way there, which 
is described by a nonequilibrium ionization model. In Xspec, the appropriate models are 



vequil for equilibrium conditions, and vnei for nonequilibrium conditions (Borkowski et al. 



2001, 1994 Hamilton et al. 1983 Liedahl et al. 1995). These models are convolved with a 



photoelectic absorption model, phabs, which accounts for interstellar absorption along the 



line of sight. Dielectronic recombination rates are taken from Mazzotta et al. (1998) with 



abundances from Anders & Grevesse (1989) and cross sections from Baluncinska-Church & 



McCammon (1992). The line list used is an improved version of the APED line list (Smith 



et al. 2001), which has been developed by Kazik Borkowski and includes more inner shell 



processes especially for Fe-L lines. Additional fitting was performed using a plane parallel 



shock model, vpshock (Borkowski et al. 2001 Sedov 1993), but this never resulted in better 



fit statistics in comparison to the vequil or vnei models. 

Investigation of emission line lists for plasmas with kT ~0.1-0.2 keV show that the 
emission is dominated by high ionization states of oxygen, nitrogen, and carbon. At temper- 
atures of kT > 0.2 keV iron and neon have significant contributions, especially above photon 
energies of 0.7 keV. The fits presented here typically varied the abundances for some sub- 
section of these five elements. All non-varying elemental abundances are set to solar values. 
Allowing Nh to vary causes it to trend toward unphysically low values with large errors in 
a handful of regions. Therefore, for all fits the column density of intervening material, Nh, 
was set to 8 X 10^° cm~^ consistent with other Cygnus Loop studies (McEntaffer & Cash 



2008 Zhou et al. 2010; Miyata & Tsunemi 2001). Other variables include the tempera- 



ture of the plasma, kT, a normalization parameter which includes emission measure, norm, 
and the ionization timescale, r. The ionization timescale is the product of the time since 
shock heating and the electron density of the shocked material, r = n^t, with > 10^^ cm~'^ s 
typically signifying equilibrium conditions. This parameter is only used in the vnei models. 



A summary of the best fit parameters is shown in Tables [T| |2} and |3j Quoted errors 
are 90% confidence intervals for the parameter of interest. The A regions are best fit with 
two temperature components, however Al requires two equilibrium plasmas while A2 uses 
two plasmas not yet in coUisional ionization equilibrium. Single temperature vequil and vnei 
fits were attempted for Al and A2 but are unable to produce a xled < 2. Substituting one of 
the equilibrium plasmas with a nonequilibrium plasma does not improve the statistic for 
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Al and the r parameter is driven to its maximum, thus suggesting equihbrium conditions. 
Region A2 could only be adequately fit {xled < 2) with nonequihbrium plasmas in its two 
component model. These fits tie the abundances between C, N, and O together for each 
component. In addition to a relatively high kT component, Al shows the presence of a low 
temperature plasma at kT = 0.09 keV while A2 is characterized by only high temperature 
gas. The abundances in all components are reduced from solar. 

In the B regions, all extracted spectra are best fit with a two temperature equilibrium 
model. It is only necessary to vary the abundance of O to obtain an adequate fit. The 
abundance of O is tied between the two temperature components. Allowing C and N to vary 
does not improve the fit statistics and if allowed to vary their values are typically driven to 
zero or to unphysically high values, e.g. tens to hundreds. Again, several combinations of 
models were attempted to fit each B region. Single temperature fits {vpshock, vequil, and 
vnei) were incapable of fitting the data from any of the regions. In region Bl only a two 
component vequil model is reasonable. Changing one of the components to nonequihbrium 
conditions does not improve statistics and the r parameter pegs high. Including a nonequi- 
hbrium component in the B2 fit does produce a reasonable fit {Xred — 1-3); but an F-test just 
favors the two component equilibrium model (F-test value = 3.1 with probability, P = 0.06, 
where Pcru = 0.05). B3 is similar to Bl in that the inclusion of nonequihbrium conditions 
does not improve the fit and the r parameter reaches its upper limit. In addition to the two 
component equilibrium model, the B4 spectrum can be fit by both a vequil + vnei model and 
a vnei + vnei model {Xred = 1-6 for both), but an F-test clearly favors the double equilibrium 
model (probabilities of 0.34 and 0.45, respectively). Finally, a vequil + vnei model fits B5, 
but with a marginally improved Xred leading to an F-test probability of 0.94, thus negating 
the addition of the nonequihbrium model. Again, the O abundance is depleted from solar 
in all regions with a range of 0.51-0.82 (where an abundance of 1.0 would signify solar). In 
each region there is a low temperature, soft component with kT ranging from 0.09-0.12 keV, 
and a harder component with temperatures ranging from 0.26-0.6 keV. 

Table |2] contains the fit parameters for the bright regions encompassed by CI and C2. 
In both cases, the given model is the only reasonable fit obtained for these regions. The C2 
fit is very similar to the B regions - O is the only varied abundance in a two component 
equilibrium model consisting of a low T and high T plasma {kT = 0.11 and 0.37 keV). The 
CI fit is similar to the A2 fit, but at lower temperatures-an average of 0.19 keV versus 0.27 
keV. The abundances of C, N, and O were allowed to vary independently, but each was tied 
between the two temperature components. The abundances of N and O are again depleted 
but there is a suspicious upper limit to C. Lines from C VI lie below our energy range and 
would contribute to the shape and trend of the extracted spectrum at low energies. This high 
C value is a result of a single high value in the lowest energy bin. Maintaining a constant 
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level of interstellar absorption leads to an increase in the C abundance to compensate for 
this high bin. 

The E regions sample the soft X-ray enhancement that extends to the west. The best 
fit parameters for these regions are given in Table [3j In contrast to the other regions, it is 
necessary to vary the Ne and Fe abundances in E1-E3 to obtain an adequate fit. This leads 
to depletions in O and Ne (0.12-0.25 and 0.32-0.41, respectively) with Fe close to or slightly 
higher than solar (0.8-1.6). Region El is adequately fit with vequil alone as well as a two 
component vequil + vnei model, but F-tests performed between these models and vnei favor 
the single temperature nonequilibrium model. E3 is much different in that only a single 
temperature equilibrium fit is valid. Attempts at using vnei, alone or in combination, drives 
the ionization parameter to equilibrium values. Finally, the only model capable of fitting 
region E2 is a two temperature component vequil + vnei model. The high temperature 
nonequilibrium component is similar to regions El and E3 in temperature and elemental 
abundances, but the additional equilibrium component is at a lower temperature with solar 
abundances. 

The extracted spectra from each of the regions discussed above are given in figure [5j In 
addition to the extracted data, the best fit model is shown for each region as a solid line. 
Regions B2 and B3 have been chosen as representative spectra from the B regions. Visual 
inspection of figure [5] shows that the B and C regions are dominated by soft emission below 
700 eV. More specifically, other than the strong O VII triplet from 561-574 eV, no significant 
emission lines stand out. Alternatively, the emission from the A and E regions appears to 
be harder or at least contains a higher temperature component. First, in addition to the 
bright O VII triplet, the presence of O VIII at 654 eV, 775 eV, and 817 eV can be inferred. 
This may solely be due to the improved count statistics in these bright regions. More 
convincingly, there are humps of significant emission above 700 eV present. This emission is 
most likely from Fe XVII lines which are expected around 725 eV, 727 eV, 739 eV, 812 eV, 
and 826 eV, as well as the He-like Ne IX triplet around 905-922 eV. The peaks of emission 
lines are systematically underestimated which may point to an error or conservatism in the 
redistribution matrices. 

Of final note, the bright knot in the southeast stands out even in this highly complex 
field. In an attempt to exploit the high spatial resolution of Chandra, we extract 42 spectra 
from this region using a fine grid of boxes, each measuring 27 x 22 pixels (Figure [6]). Given 
a distance of 540 pc, each box probes a mere 0.034 x 0.028 parsecs of area at the remnant. 
The extracted spectra were fit with the same round of models but only the two temperature 
component models provide an adequate fit. In general, these regions are well fit with a 
vequil + vnei model. Of the 42 regions, 7 were fit statistically better (via F-test) with a 



- 8 - 



double nonequilibrium model. Parameter maps of the best fit parameters for kT and rie over 
the grid are shown in figure [7] where is calculated from the fit normalization. The lower 
temperature component corresponds to kTi and Uei and is typically in equilibrium while the 
higher temperature component has kT2 and ne2 and is always in nonequilibrium. The average 
error in the equilibrium temperature is ~ 3% and for the nonequilibrium temperature it is 
~ 4%. The average errors in density range from 8-11%. The abundances of C, N, and 
O were allowed to vary but fixed to each other. The abundances were typically depleted 
relative to solar to a similar degree as the other regions (0.3-0.8 times solar) with a few grid 
boxes around solar values. The xled average and standard deviation of 1.5 ± 0.4 over 

the grid. 



4. Interpretation 

In general, the extracted spectra are best fit by two temperature component models. 
A handful of scenarios may be used to explain this phenomenon. First, the region may be 
merely confused due to multiple, physically distinct plasmas along the line of sight. This 
effect is expected to become less of an issue toward the rim and more significant for the 
interior regions. Second, the region may encompass a high density region and a low density 
region, both of which have been recently shocked. The high density region will equilibrate 
(larger ionization parameter, riet) and cool more rapidly than the low density region, thus 
leaving two temperature components. Finally, a similar situation may exist with two distinct 
density regions, but regions that are causally linked. Such a situation may occur when the 
shock wave passes out of the precursor formed cavity and into the cavity wall. The high 
density gas in the cavity wall will again equilibrate more rapidly, but the low density region 
is located interior to the wall and a reverse shock will reshock and additionally heat this 
interior gas. 

In order to distinguish between these scenarios (or develop a new one) we calculate 
the densities (rze) for each of the model components in all regions. The results are given 
in Table |4| The density is calculated using the norm parameter within each component 
where norm = [10~^'^/{4ttD'^)] J nenndV. A distance of 540 pc is used for D and it is 
assumed that n^ = nu- The volume is calculated using the area of each region multiplied 
by the depth of the column of remnant gas each region intersects. This depth is found by 
assuming the remnant is a filled sphere with radius 1.5° or 14 pc and calculating a chord 
length based on the radial position of each region. It is assumed that the B regions sample 
the outermost radii of the sphere. These assumptions probably lead to the error in volume 
being the largest contributor to the error in density so that the densities in table |4] should 
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be treated as approximate. However, since the volume between components within a region 
remains constant, the relative density between components in a given region is accurate. 

The most evident trend is in the B regions. In all 5 regions the lower temperature 
component {{kT) ~ 0.10 keV) exhibits a significant density enhancement - on average, over 
an order of magnitude higher than the high temperature component {{kT) ~ 0.38 keV). This 
result, combined with the location of the B regions at the easternmost edge of the X-ray 
emission strongly suggests that these B regions trace out the location of the cavity wall. The 
high density, low temperature component represents the high density enhancements in the 
cavity wall itself while the low density, high temperature gas most likely lies just interior 
to the wall in the precursor formed cavity. The blast wave has shock heated high density 
clouds in this wall and reverse shocks from this interaction propagate toward the interior. 
The high density gas equilibrates rapidly, cools via emission lines (predominantly O VII) and 
has a lower temperature than the lower density interior which has been reshocked and takes 
longer to cool. 

It may be expected that the interior lower density gas is in nonequilibrium, but that 
is only true if it has been very recently shocked, i.e. t and Ue are sufficiently low to give a 
low ionization parameter. However, in all B regions the high temperature component can 
never be fit with a nonequilibrium plasma and any attempts drive the ionization parameter 
to its upper limit, which signifies equilibrium and therefore a longer timescalc since being 
shocked. Solving the Rankine-Hugoniot relations for a strong shock and 7 = 5/3 gives the 
post-shock temperature as a function of shock velocity, T = {3/16)fiv'^/k with /i = 1.2. This 
can be easily applied to the cavity wall where the average shock velocity in this high density 
gas is 209 km/sec. If we then assume a Sedov solution the time since shock heating can 
be calculated as t — {2/5)R/v where R is the radial shock location, 14 pc. While this is 
an inaccurate assumption in the cavity wall it will provide an upper limit to the time since 
the shock velocity in the dense gas is lower than the Sedov velocity. This upper limit is 
~26,000 years. Ideally we could use the shock velocity combined with visible spatial offsets 
between the shocks in the cloud and the reverse shock to estimate a propagation timescale, 
but the resolution and/or count rate is not present in this observation. The situation for 
the interior gas is even more complicated if we assume it has been shocked twice, once by 
the initial supernova blast wave and again by a reverse shock. The contribution from each 
of these shocks to the temperature of the gas is uncertain. Therefore, a similar calculation 
will overestimate the shock velocity. Even so, the high temperature component in these 
regions is still best fit with an equilibrium plasma which suggests high densities and/or 
long timescales since shock heating. Therefore, this crude age estimation, the more precise 
velocity calculation, and the large ionization timescales suggest that the B regions have 
equilibrated due to their high density and possibly also due to a significant time since being 
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shocked. 

In other regions the situation is comphcated by their interior locations and presumed 
hne of sight confusion. For instance, the A regions exhibit similar properties to the B 
regions with high densities typifying the low temperature component. This effect is more 
pronounced in region Al given the more distinct temperature difference between components 
in comparison to region A2. However, given the location of these regions it is more likely 
that line of sight confusion leads to multiple components rather than probing a single distinct 
cavity wall/shock interaction as with the B regions. Inspection of figure |4] shows that Al 
is dominated by hard emission, but may lie close enough to the rim of the remnant to pick 
up some low temperature cavity wall emission along the line of sight. This is supported by 
the presence of a very high density, low temperature component in the Al fit. Also, the 
presence of soft x-rays is somewhat evident in both the false color and ratio maps. In region 
A2 though there is a distinct absence of soft X-ray emission in figure |4| The A2 line of sight 
is located interior to Al and therefore probes through a much larger column of the high 
temperature gas, which dominates as shown by the best fit temperatures of kT ~0. 23-0. 30 
keV. If we assume that the low temperature component is completely contained within a 
shell of thickness comparable to a B region, and use a radius of 14 pc for the location of the 
cavity wall, then we can calculate the relative contributions of the high and low temperature 
components in each of these A regions. Region Al probes a total shell volume of 0.015 pc^ 
and interior volume of 0.12 pc'^ while the Al line of sight contains a shell volume of 0.011 pc'^ 
and interior volume of 0.16 pc^. This estimates that region A2 contains approximately twice 
as much interior, high temperature emission relative to low temperature emission. The lower 
density, high temperature plasmas of A2 are consistent with nonequilibrium conditions for 
this interior region. Furthermore, the interior plasmas of these A regions display significant 
emission from higher energy lines, most likely Fe XVII and Ne IX, also consistent with the 
high temperatures in these regions. Therefore, while the A regions probably encompass 
some cavity wall emission, some reverse shock emission, some singly shocked low density hot 
emission, or some combination of all of the above, it is unclear if the temperature/density 
components are related given that these regions probe an interior sightline that is prone to 
confusion. 

The C regions follow a similar trend as the A regions. C2 shares properties of Al; the 
line of sight probes both a low temperature, equilibrium plasma near the edge of emission 
and hotter ineterior gas as well. The presence of the low temperature component is quite 
evident upon examination of the ratio map. C2 outlines the brightest region in this map so 
it is not surprising that there is a significant soft component. This component is probably 
due to a density enhancement in the cavity wall (rig = 5.6), but possibly located in the 
foreground, toward the observer in comparison to the B regions, assuming a spherical shape 
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to the remnant. However, the exact location of the soft emission along the line of sight is 
complicated by the high degree of complexity in this area as seen in the raw data and color 
image. Region CI gives no further clues either given the distinctly different fit parameters 
in comparison to C2. This particular line of sight probes mostly higher temperature gas. 
It is quite bright in the raw data, but shows no significant abudance of soft X-ray emission 
in the ratio map, especially in comparison to C2. The lack of soft X-rays leads to the 
relatively higher fit temperatures {kT = 0.19 keV versus 0.10 keV on average for soft X-ray 
components) while the bright emission is consistent with the high density, i.e. large emission 
measure, of this region. Region CI is the only example of a high density, high temperature, 
nonequilibrium plasma in any of the extracted regions. If this density enhancement is in 
nonequilibrium then it must have been recently shocked. Using the ionization parameter 
and the calculated density show that this plasma was most likely shocked within the last 
1000 years. It is possible that this cloud traces out where the blast wave is interacting with 
the cavity wall but at a larger radius than the B regions, thus leading to the more recent 
shock time. If this is the case, then the cloud must be located considerably in the foreground 
relative to B. Therefore, the C regions are again confused by line of sight projection yet 
show the presence of probable cavity wall emission in addition to the hotter, lower density, 
interior emission. 

The picture is not so clear in the E regions. These extractions are even farther toward 
the interior than the A regions, but seem to have simpler fits. For instance, a single tem- 
perature equilibrium plasma is the only reasonable fit obtained for region E3. One might 
expect that a multiple component fit at high temperature and nonequilibrium may be fa- 
vored. However, these regions were chosen because they sample the soft X-ray extension 
that protrudes westward from the edge of the remnant (see Figure |4]). It may be the case 
that region E3 looks directly at a dense cavity wall cloud, which obscures the view into the 
hotter, low density interior. The calculated density of 3.0 cm~^ is not overly dense, but 
still a few-several times denser than interior densities seen in other regions. Region E2 on 
the other hand is well fit with a two component model containing both a low temperature 
equilibrium plasma with a high temperature nonequilibrium plasma. This again may be in- 
dicative of cavity wall emission in a dense medium that has equilibrated and cooled combined 
with interior plasma that is at a higher temperature (possibly reshocked) and lower density 
leading to nonequilibrium conditions. However, the calculated densities are not as clear cut 
as they are in previous regions. The low temperature component for E2 is at approximately 
the same density as the high temperature component thus lending no resolution to the line 
of sight confusion. Inspection of the color map shows that while E3 is situated in the midst 
of soft emission, E2 is sampling a region with a mix of low and higher energy emission with 
no obvious distinction. Finally, El is even further afield from the soft emission of E3 than 
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E2 and appears to be dominated by harder X-rays. However, the best fit spectrum is only 
marginally hotter than E3 at 0.17 keV and at nearly the same density. The main difference 
between these regions is that El appears to be in nonequilibrium and therefore more recently 
shocked relative to E3. Perhaps El is a medium density cloud that has been reshocked by a 
reverse shock originating in the denser cavity wall. While the E regions provide interesting 
hints into the morphology and possible interactions in this section of the Cygnus Loop, it 
is impossible to state anything definitively given the confusion inherent to these interior 
regions. 

The grid analysis of the bright southeast knot shown in Figure [7] has presented some 
intriguing results as well. Inspection of the equilibrium maps of fcTi and riei reveal that they 
appear almost as negatives of one another. As expected, the highest densities correspond to 
areas of bright emission with density falling off toward the edges of the grid. These dense 
pockets have slowed the shock considerably more than the surrounding shock regions. The 
slower shock velocity results in correspondingly lower post-shock temperature. The densities 
calculated here are similar to those found in the cavity wall along the B regions. This 
bright area of emission could be an indicator of a protrusion into the cavity wall or mark 
a particularly large, dense cloud along the cavity wall slightly in the foreground relative to 
the B regions, or a combination of the two. Outside of the bright emission the color map 
contains harder emission which is consistent with the higher temperature fits in these grid 
squares. These regions could be located outside of the dense cloud and probe interior gas in 
the background, along the line of sight. The brightest knot in the southeast, which displays 
high density, is located close to the edge of the X-ray emission. In contrast to region B 
which appears to follow an ~0.43 pc vertically elongated section of the cavity wall, this 
particular emission seems to be located in a more discrete cloud. This cloud could possibly 
be oriented along the line of sight, which may help in explaining the intensity of emission 
and the apparent indentation in the X-ray emission at this point. Previous studies of the 
XA region have explained the X-ray emission as arising from a single, large protrusion which 
is responsible for the entire field of view studied here. However, the fine spatial resolution 
of Chandra has uncovered much smaller blast wave/cloud interactions, possibly as small as 
a single grid square, 0.035 pc x 0.028 pc. 



5. Conclusions 

The high spatial resolution of Chandra has given us insight into the structure of the 
complicated XA region. Specifically, the chain of B regions appear to follow the location of 
the blast wave interaction with the cavity wall. These B regions are located at the very edge 
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of X-ray emission and are dominated by soft X-ray emission. Spectral analysis has shown 
low temperature gas at densities an order of magnitude higher than more recently shocked 
high temperature gas. Shock heated gas in the high density cavity wall has equilibrated 
rapidly and is efficiently cooling through O VII emission. A reverse shock has propagated 
from the initial supernova blast wave interaction with the cavity wall and has reheated the 
interior gas. Furthermore, the time since interacton of the reverse shock with the interior 
gas has been long enough for the plasma to equilibrate as well. Given this interpretation, 
the precursor wind blew a cavity that is ~ 14 pc in radius. 

The identification of the cavity wall is simpler and more robust than identifying struc- 
ture interior to the edge. While the color map shows a general trend of soft X-ray emission 
on the exterior with harder emission toward the interior, this trend is not obvious in discrete 
extractions. Regions interior to the B regions have line of sight confusion and significant 
chord length through the spherically shaped remnant. Some of these regions show strong 
signs of including the cavity wall (high density, low temperature components), but the in- 
terpretation is less obvious given that these regions are not located at the rim of emission. 
Furthermore there is a mix of temperatures, densities, and ionization states with no strong 
trends. However, we feel that this is a significant finding given that previous studies support 
a much simpler morphology whether it be a single finger-like protrusion into this region or a 
few clumps. We would argue that the fine spatial resolution here has uncovered a situation 
where any region removed from the edge of emission will have significant line of sight con- 
fusion and may or may not include cavity wall emission, interior emission, interior reverse 
shock emission, or multiple components of each. It is not to say that any given region is 
physically disconnected from another since it is definitely the case that the X-ray morphology 
is arising from a global interaction of the blast wave with a precursor formed cavity wall. 
However, the intricate details of the X-ray emission and therefore the physical parameters of 
the plasma can vary over the smallest scales. This effect is evident in the gridded extraction 
region around the bright knot of emission in the southeast. The calculated densities in the 
grid show the presence of an isolated blast wave/cloud interaction possibly dominated by 
only a few grid sections. These data display a significant increase in resolution over the 



previously highest resolution study (Zhou et al. 2010) where extraction regions are larger 
than the entire grid. 

The Cygnus Loop is clearly complex both spatially and spectrally. Chandra can image 
fine structures indicative of distinct interaction regions, but the nearly broadband spectra are 
lagging in quality. The next step in understanding the intricacies inherent to the interaction 
of this supernova with the ISM is to obtain high quality spectral data. Higher spectral 
resolution is required to constrain the parameters of the models and test their validity. 
Grating instruments onboard existing X-ray telescopes are capable of high spectral resolution 
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but are designed for point sources. The diffuse nature of galactic supernova remnants cause 
current grating observations to be signal limited with blurred spectral resolution. To increase 
understanding of these objects future instruments will need to be capable of performing high 
efficiency and high spectral resolution observations of diffuse soft X-ray sources. 

6. Acknowledgments 

The authors would like to acknowledge internal funding initiatives at the University of 
Iowa for support of this work. The data used here were obtained from the Chandra Data 
Archive. 

REFERENCES 

Anders, E., & Grevesse, N. 1989, Geochim. Cosmochim. Acta, 53, 197 

Arnaud, K. A. 1996, in ASP Conf. Ser. 101, Astronomical Data Analysis Software and 
Systems V, ed. G. H. Jacoby, & J. Barnes, 17 

Balucinska-Church, M., & McCammon, D. 1992, ApJ, 400, 699 

Blair, W. P., Sankrit, R., & Raymond, J. C. 2005, AJ, 129, 2268 

Borkowski, K. J., Sarazin, C. L., & Blondin, J. M. 1994, ApJ, 429, 710 

Borkowski, K. J., Lyerly, W. J., & Reynolds, S. P. 2001, ApJ, 548, 820 

Danforth, C. W., Blair, W. P., Raymond, J. C. 2001, AJ, 122, 938 

Fesen, R. A., Blair, W. P., Kirshner, R. P. 1982, ApJ, 262, 171 

Fruscione, A., et al. 2006, Proc. SPIE, 6270, 62701V 

Hamilton, A. J. S., Chevalier, R. A., & Sarazin, C. L. 1983, ApJS, 51, 115 
Hester, J. J., & Cox, D. P. 1986, ApJ, 300, 675 

Katsuda, S., Tsunemi, H., Kimura, M., & Mori, K. 2008, ApJ, 680, 1190 
Leahy, D. A. 2004, MNRAS, 351, 385 
Levenson, N. A. et al. 1997, ApJ, 484, 304 



- 15 - 

Levenson, N. A., Graham, J. R., Keller, L. D., Richter, M. J. 1998, ApJS, 118, 541s 

Levenson, N. A., Graham, J. R., Walters, J. L. 2002, ApJ, 576, 798 

Levenson, N. A., & Graham, J. R. 2005, ApJ, 622, 366 

Liedahl, D. A., Osterheld, A. L., Goldstein, W. H. 1995, ApJ, 438, L115 

Mazzotta, P., MazziteUi, G., Colafrancesco, S., & Vittorio, N. 1998, A&AS, 133, 403 

McCray, R., & Snow, Jr., T. R 1979, ARA&A, 17, 213 

McEntaffer, R. L., & Cash, W. 2008, ApJ, 680, 328 

Miyata, E., & Tsunemi, H. 2001, ApJ, 552, 624 

Miyata, E., Katsuda, S., Tsunemi, H., Hughes, J. P., Kokubun, M., & Porter, F. S. 2007, 
PASJ, 59, 163 

Sankrit, R., Blair, W. P., Cheng, J. Y., Raymond, J. C, Gaetz, T. J., & Szentgyorgyi, A. 
2007, AJ, 133, 1383 

Sedov, L. 1993, Similarity and Dimensional Methods in Mechanics (10th ed.; Boca Raton, 
Fla. : CRC Press) 

Smith, R. K., Brickhouse, N. S., Liedahl, D. A., Raymond, J. C. 2001, ApJ, 556, L91 

Spitzer, L. 1998, Physical Processes in the Interstellar Medium (Wiley Classics Library 
Edition; New York, NY; Wiley Interscience) 

Zhou, X., Bocchino, P., Miceh, M., Orlando, S., & Chen, Y. 2010, MNRAS, 406, 223 



This preprint was prepared with the A AS IM^jX macros v5.2. 



- 16 - 



Table 1: Best Fit Parameters for Regions A and B. 



Regioij^ 




CNOj^] 




nornii 


m 


CNO2 


1"2 


norni2 








(keV) 




(10"crn"'\s) 


(io-*/(] 


(keV) 




(10"cm->\s) 








Al 


0.086llS:S 


0.109±0.006 


n/a 


280±10 


0.236±0.003 


0.69±0.04 


n/a 


2.54±0.08 


7100 


49.4/35 


A2 


r, 909+0.005 

U.ZOZ_Q QQg 


0.6±0.1 


n 0033+"""^ 

U.UUOO_Q Q0Q3 


29±4 


0.299±0.003 


0.61±0.02 


0.60±0.05 


1.94±0.04 


9391 


51.6/40 


Bl 


0.095±0.001 


0.59±0.03 


n/a 


12.8±0.7 


n 07+0.05 


0.59±0.03 


n/a 


0.07±0.01 


1275 


24.1/13 


B2 


0.085±0.001 


0.51±0.03 


n/a 


25±1 


0.26±0.02 


0.51±0.03 


n/a 


0.17±0.03 


1176 


21.6/13 


B3 


0.096±0.001 


0.82±0.06 


n/a 


5.9±0.4 


0.27±0.02 


0.82±0.06 


n/a 


0.12±0.02 


961 


21.7/13 


B4 


0.116±0.002 


0.55±0.04 


n/a 


4.3±0.3 


n 00+" '"' 

U-5O-0.08 


0.55±0.04 


n/a 


0.03±0.01 


951 


17.1/11 


B5 


0.123±0.002 


0.79±0.05 


n/a 


2.7±0.2 


().6±0.1 


0.79±0.05 


n/a 


0.037±0.009 


1061 


26.2/14 



Note. — Abundances are shown relative to solar with solar values equal to 1. 

"All B regions vary only O. C and N are set to solar values. CNOi and CNO2 are 
independent m A regions but fixed equal in B regions. 

^Total counts per extraction region 

'^Normalization parameter, where A = [10~"^'V (47rZ)^)] j n^nndV, D is the distance to the 
Cygnus Loop and the integral is the emission measure. 
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Table 2: Best Fit Parameters for Region C. 



Regioi|] 


(keV) 


Ti nornii 
(10"cm-''s) (lO-*.-^ 


kT2 
(keV) 


T"2 

(10"cm-'ls) 


11011112 

(10-^.4) 


C 


N 







xVdof 


CI 


0.179±0.006 


0.06l|J;UJ 0.66±0.07 


0.191±0.002 


0.0043±0.0007 


38±2 


<7.5 


0.7±0.2 


0.81±0.05 


2576 


27.4/16 


C2 


0.107±0.001 


n/a 11.7±0.5 


S7+" °'^ 


n/a 


0.04±0.01 


1.0 


1.0 


0.51±0.03 


1597 


25.8/18 



Note. — Abundances are shown relative to solar with solar values equal to 1. 
"Total counts per extraction region 

^Normalization parameter, where A = [10""^'*/ (AttD'^)] J nf.riHdV , D is the distance to the 
Cygnus Loop and the integral is the emission measure. 
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Table 3: Best Fit Parameters for Region E. 



Regioij] 


kTi 


Ne 


Fe 


n 


nornii 




norm2 




xVdof 




(keV) 






{10"cm-''.s) 


(10-"/^ 


(keV) 


{ur"A) 






El 


0.174±0.001 


0.217±0.009 0.35±0.05 


0.8±0.1 


o r+0.5 
^■■-•-0.4 


10.5±0.3 


n/a 


n/a 


4333 


50.7/28 


E2 


0.198±0.002 


0.12±0.01 0.32±0.07 


1 


0.8±0.1 


4.6±0.2 


0.146±0.003 


2.4±0.1 


3910 


55.3/27 


E3 


0.166±0.001 


0.25±0.01 0.41±0.07 


1.6±0.2 


n/a 


12.1±0.3 


n/a 


n/a 


4330 


45.8/31 


Note. 


— Abundances 


are shown relative to solar 


with sol 


ar values equal to 1. 











"Total counts per extraction region 

^Normalization parameter, where A = [10~-^'"Y (47rD^)] J neTijjdV, D is the distance to tlie 
Cygnns Loop and the integral is the emission measure. 
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Table 4: Electron Density and Shock Velocity. 



Regior 


rie (cm ^) 


Vol. (pc^) 


V (km/sec) 


Al 


13, 1.2 


0.041 


192, 317 


A2 


3.4, 0.8 


0.058 


315, 357 


Bl 


9.5, 0.7 


0.012 


201, 397 


B2 


13, 1.0 


0.012 


190, 332 


B3 


6.4, 1.0 


0.012 


202, 339 


B4 


5.4, 0.4 


0.012 


222, 403 


B5 


4.4, 0.6 


0.012 


229, 506 


CI 


1.4, 10 


0.010 


276, 285 


C2 


5.6, 0.4 


0.010 


272 


El 


2.8 


0.035 


250, 291 


E2 


1.4, 1.8 


0.035 


192, 317 


E3 


3.0 


0.035 


266 



Note. — The first density listed corresponds to the lower T component for regions with 
two components. 
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Fig. 1. — Image of the raw data from the Chandra ACIS S3 back-illuminated CCD. The 
pixels have been binned by 2 and are shown linearly scaled with values ranging from 0-70 
counts per bin. 



35.0 30.0 25.0 20.0 15.0 10.0 05.0 20:57:00.0 
Right ascension 



Fig. 2. — A false color image of the raw data. The three colors, red, green, and blue, depict 
the total counts in three energy bands, 0.3-0.5 keV, 0.5-0.8 keV, and 0.8-1.5 keV, respectively. 
The direction of shock propagation is from right to left so that the eastern red edge is closest 
to the current shock position. 
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35.0 30.0 25.0 20.0 15.0 10.0 05.0 20:57:00.0 
Right ascension 

Fig. 3. — The ratio of detected soft energy flux (0.3-0.5 keV) to hard energy flux(0.8-1.5 
keV). The bright regions trace the soft X-ray emission arising from the blast wave encoun- 
tering the local ISM. 
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Fig. 4. — Spectral extraction regions are overlaid on the band ratio map and the false color 
image. The regions are labeled are on the right and were chosen based on the morphology 
evident in these images. 




Fig. 5. — Spectra from the extracted regions with best fit model shown as sohd hnes. The A 
regions are best fit with a nonequihbrium model, the B and C regions with an equilibrium 
model, and the E regions with a two temperature component model. 



Fig. 6. — Close-up of the bright knot of emission in the southeast of the S3 CCD field of 
view. The bottom image shows an overlay of the extraction region grid. 
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kT, (keV) kT; (keV; 




Fig. 7. — Best fit parameter maps for fits to the bright knot region in the southeast of this 
observation. A spectrum was extracted and fit for each of the 42 boxes shown which corre- 
spond to the grid of boxes in figure |6] Fit parameters include temperature and density for 
each of the two model components. The lower temperature, typically equilibrium component 
has parameters kTi and n^i with kT2 and characterizing the nonequilibrium component. 
The vertical color bar gives the corresponding range of values for the colors in the boxes. 



